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Abstract 

Superluminal neutrinos are expected to lose energy due to bremsstrahlung. It 
is dominated by e + e~-pair production if kinematically allowed. The same signature 
was used in searches for 3-body decays of hypothetical heavy sterile neutrinos. From 
the published analyses of these searches performed by CERN PS191 and CHARM 
experiments we set upper limits on the neutrino velocity in the energy range from 
0.2 GeV to 280 GeV. Our limits are well below the neutrino velocity favored by the 
recent OPERA results. For energy-independent neutrino velocity the limits obtained in 
this paper are stronger than those coming from ICARUS experiment and observations 
of Supernova SN1987a. 

Recently neutrino physics has attracted renewed attention because of the surprising result 
of the OPERA collaboration [1] (see, however [2]). Their claim is that the travel time of 
muon neutrino with average energy 17.5 GeV is smaller than the light travel time at the 
confidence level above 5 a. This result was interpreted as a superluminal propagation of 
muon neutrino, which can be parametrized by 



where v v is the neutrino velocity and we set the Planck constant H and photon velocity c 
equal to unity, h = c = 1. 

However, Cohen and Glashow [3] put forward a serious objection to this superluminal 
interpretation. They considered several decay processes which would take place in vacuum 
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5 = (V - 1) ~ 5 x icr 5 , 




under the assumption 5 > 0. It was pointed out that by analogy to the Cherenkov radiation, 
the electroweak process 

^ ->■ + e + + e~ (2) 

should occur, and the neutrino spectrum should be drastically changed as neutrinos travel a 
distance of about 730 km from the target at CERN to the OPERA detector at Gran Sasso. 
The neutrino energy spectrum near the OPERA detector was measured by the ICARUS 
collaboration, and no attenuation was observed. Using the result of Cohen and Glashow [3] 
and the analysis of the ICARUS data, the limit 

5 < 4 x 1(T 8 (3) 

was obtained [4] under the assumption of energy-independent 5. 

In this paper we give interpretation of other neutrino experiments, where the emission 
of e + e~-pair could be observed, in terms of the decay (2). The rate of the decay (2) in the 
laboratory frame is given by 

r = V^aW , (4) 

where Gp is the Fermi constant, E is the neutrino energy and k' is a numerical constant. 
In the original paper [3] the value k' = 1/14 was used for estimates. As shown in Ref. [5], 
this value is model-dependent and can be made two or three times smaller. Due to the very 
strong dependence of (4) on 5 and neutrino energy the precise value of k' is not crucial for 
our estimates and we adopt k! = 1/14, which was also accepted in Ref. [4] when obtaining 
the upper bound (3). The threshold of the pair production (2) is 

E = 2m e /VS . (5) 

The signature of the process (2) is the appearance of e + e~-pair from nothing downstream 
of a neutrino beam. Remarkably, the same signature is exploited in searches for heavy 
sterile neutrino decays. Indeed, sterile neutrinos mix with the Standard Model (SM) active 
neutrinos thereby producing masses and mixing in active neutrino sector required to explain 
the neutrino oscillations. This mixing also gives rise to sterile neutrino production in meson 
decays and subsequent sterile neutrino decays into SM particles including e + e~z/, i.e. the 
same final state as in (2) (for a concrete model example see, e.g., Ref. [6]). 

In the neutrino energy range E ~ 1-8 GeV the stringent bound on the lepton pair pro- 
duction follows from the analysis of CERN PS191 experiment. Neutrinos were eventually 
produced on a beryllium target by protons of energy 19 GeV. This experiment was specially 
designed [7] to look for possible sterile neutrino decays to e + e~u in vacuum and has observed 
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no events with collected statistics of 8.6 x 10 18 protons on target. The differential energy 
spectrum of muon neutrinos dN u /dE expected from Monte Carlo simulations in the range 
of energies from 200 MeV to 8 GeV is presented in Ref. [8] and shown in Fig. 1. In the case 
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Figure 1: Neutrino spectra in CERN PS 191 and CHARM experiments. 



of process (2), neutrino loses about three-quarters of its energy in each emission [3], so even 
for E ~ 200 MeV the energy of leptons is well above the ionization threshold of the detector, 
and the efficiency of registration was close to 100% [8]. 

For the vacuum camera of length / w 10 m each neutrino has a probability of e + e~-pair 
production 

P = P (£ ) = i.r«o.8xio-"x( T A s ) 3 ( T1 | v ) 5 . (6) 

For a given neutrino energy interval one writes for the number of e + e~-pairs from nothing: 

N = JdE d ^P(E). (7) 

In the original search for sterile neutrino [7] the obtained limit was based on a set of require- 
ments on the signal events, one of which is the absence of any double track events due to the 
e + e~ pair. In our case the angle between electron and positron is small, ~ y/~5, well below 
the detector angular resolution, ~ 10~ 2 , see Fig. 2. To constrain neutrino superluminality we 
use the subset DST2 of ~ 4000 events which (judging on description presented in Ref. [9]) 
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includes our signal events (if any): one expects a single track with double heat in electro- 
magnetic calorimeter. Having no real data in hand, we cannot proceed further in analysis 
of DST2. Then, to obtain the limit on 5, we do not subtract SM background and suppose 
that all these ~ 4000 are due to the Cohen-Glashow effect. One can set an upper limit on 
5 at 95% CL. by adopting the normal distribution and requiring N < 4100. Then for the 
energy-independent 5 we integrate in r.h.s. of Eq. (7) over the entire interval of neutrino 
energies 0.2-8 GeV and obtain 

5 < 2.2 x 10" 6 at 95% CL . 

In the case of energy-dependent 5 = 5 (E) one has to integrate in r.h.s. of Eq. (7) to obtain 
the limits on neutrino dispersion relation. In particular, if 5 (E) does not vary significantly 
in the energy interval E ± E/2, one can exploit the relation 

r E+E/2 Jju n2 pE+E/2 inr 

4100 > N E = / dE' d -^ P(E') ~ 15* (E) h>±^ / dE> ^ E' 5 
Je-e/2 d E 192 7r d J e _ e/2 dE> 

to set a limit on moderately-varying 5 (E) as presented in Fig. 2. 
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Figure 2: Upper limits on 5 (E) from the results of CERN PS191 and CHARM experiments; 
straight line refers to the pair-production threshold (5), so the region below this line is 
kinematically forbidden. 
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We obtain similar limits from the negative results of CHARM experiment [10] which also 
performed searches for sterile neutrino decays. In this experiment, the CHARM detector of 
length I fa 12 m was exposed to wide-band neutrino beam produced by ~ 1.4 x 10 18 protons 
with energies ~ 400 GeV incident on a copper target. The neutrino spectrum expected from 
Monte Carlo simulations 1 is given in Ref. [11] and shown in Fig. 1. The selected events 
were required to have a shower energy W deposited in the CHARM calorimeter between 7.5 
and 50 GeV and a value of the variable W 2 8 2 below 0.54 GeV 2 (6 is the angle between the 
shower axis and the direction of incoming neutrino). A total of 331 events were selected. 
The number of events attributed to heavy neutrino decay is compatible with zero, 1 ± 41 
event. In the case of decay (2) one has 8 2 ~ 5, so the selected above cuts do not interfere 
in our analysis. To place a limit on 5 in this case we require to have less than 85 events 
due to the bremsstrahlung (2). Also, the statistics for e + e~ production from nothing was 
approximately doubled when detector was exposed to muon antineutrino flux produced by 
~ 5.7 x 10 18 protons on target [10]. For energy-independent 5 we obtain from (7) after 
integrating over the whole available neutrino energy range 10-280 GeV, and taking into 
account the 44% efficiency [10] of e + e~-registration, 

5 < 3.6 x 10" 9 . 

This limit is stronger than the constraint (3) obtained by the ICARUS collaboration. Strik- 
ingly, this limit is even comparable to the stringent direct bound [12] obtained from the 
observations [13] of neutrino signal from Supernova 1987a, 

5 < 4 • 10~ 9 , 

which is valid for the energy range around 10 MeV. However, our limit is weaker than what 
is expected [3] from the analysis of 10-100 TeV neutrino events measured by the IceCube 
experiment. The reason is the strong dependence of the emission rate (4) on neutrino energy. 

For a moderately-varying S (E) we follow the procedure used for PS191 and obtain limits 
presented in Fig. 2. 

To summarize, in this paper we have interpreted the negative results of searches for 
sterile neutrino decays in CERN PS191 and CHARM experiments as the restriction on e + e~- 
pair production by possibly superluminal muon neutrinos. We have set the upper limits on 
neutrino velocity in the energy ranges 0.2-8 GeV and 10-280 GeV supposing neutrino velocity 

1 We use wide-band neutrino beam spectrum presented in [11] for CDHS detector, which was placed right 
in front of the CHARM detector along the neutrino beam. The effective areas of these detectors were almost 
identical. We thank A. Rozanov for this suggestion. 
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to be only moderately varying with energy. In a particular model the limit can be refined by 
convolution of neutrino production rate with neutrino beam spectra according to formulas 
presented in the paper. Covering the gap between 8 and 10 GeV should be possible with 
information on high energy end of the neutrino spectrum in CERN PS191 experiment. In 
both experiments further strengthening of the bound might be possible with adjustment of 
the selected cuts. In particular with lower cut on W 2 6 2 in CHARM experiment. Potentially 
interesting limits one can expect from reanalysis of data collected at 35-meter length detector 
with CHARM calorimeter module installed, which was designed to search for heavy neutrinos 
produced by charmed meson decays [10]. 

Analogous limits are expected to be obtained by the dedicated analysis of the data from 
NOMAD experiment, see Ref. [14] for the preliminary bound. 

We thank V. Rubakov and S. Sibiryakov for discussions, and A. Rozanov and F. Vannucci 
for correspondence. 
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